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Existing coherent network analysis techniques for detecting gravitational-wave bursts simulta- 
neously test data from multiple observatories for consistency with the expected properties of the 
signals. These techniques assume the output of the detector network to be the sum of a station- 
ary Gaussian noise process and a gravitational- wave signal, and they may fail in the presence of 
transient non-stationarities, which are common in real detectors. In order to address this problem 
we introduce a consistency test that is robust against noise non-stationarities and allows one to 
distinguish between gravitational- wave bursts and noise transients. This technique does not require 
any a priori knowledge of the putative burst waveform. 
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I. INTRODUCTION 

Gravitational-wave bursts (GWBs) are among the 
most exciting classes of signals that large-scale, broad- 
band interferometric gravitational-wave observatories 
will attempt to detect. These instruments have already 
started to collect data, and are beginning to coordi- 
nate their activities in order to perform joint, world- 
wide, searches. The United States LIGO [1] observatory 
has reached its design sensitivity goal, and is currently 
conducting its fifth data-taking run. Both the British- 
German GEO600 [2] interferometer and the Japanese 
TAMA [3] detector have conducted multiple data-taking 
runs, often in coincidence with LIGO. The French-Italian 
VIRGO [ ] detector is approaching its operational phase 
with an anticipated sensitivity comparable to that of the 
US detectors. 

Potential sources of GWBs include merging compact 
objects [5, 6, 7], core-collapse supernovae [8, 9, 10, 11], 
and gamma-ray burst engines [12]. By considering the 
anticipated strengths of the gravitational-wave bursts 
emitted by these systems and the present performance of 
the detectors, it is clear that coincident experiments are 
necessary to maximize the chances of successfully identi- 
fying an astronomical event. This is because signals emit- 
ted by these sources cannot currently be modeled with 
sufficient accuracy to distinguish them from transient 
noise non-stationarities ( "glitches" ) affecting the data of 
the detectors. In addition, a real burst signal may not be 
much larger than the noise levels in the detectors. More 
detectors at the same site provide one way for increasing 
the signal-to- noise ratio (SNR), but a network of detec- 
tors spaced across the globe inherently provides much 
more information about the wave and its source. With 



a network of three or more interferometers, for instance, 
one has in principle sufficient information for good SNRs 
to infer the direction to the source. If the events are not 
corroborated by an electromagnetic detection (optical. 
X-ray, or radio) this could be crucial. Such a network 
also has enough information to reconstruct the two in- 
dependent polarization amplitudes of the wave. This is 
possible because three interferometers provide three in- 
dependent measures of the gravitational wave (which are 
functions of time) and two independent time delays. If we 
consider the geometrical plane passing through the three 
sites of the interferometers, the two independent time de- 
lays jointly identify two possible points in the sky where 
the signal could have come from that are mirror-images 
of each other with respect to this plane; see Figure 1. 
Since the detector antenna patterns are not symmetric 
with respect to this plane, it is further possible to resolve 
this two-fold ambiguity by properly accounting for the 
antenna pattern asymmetry in the analysis of the data. 
Once the sky position has been determined it is straight- 
forward to make a minimum- variance estimate of the two 
polarization waveforms of the GWB as linear combina- 
tions of the detector data streams. The extraction of this 
information from the combined responses of the individ- 
ual detectors of the network is called the "solution of the 
inverse problem" in gravitational-wave astronomy. 

The inverse problem for gravitational- wave bursts with 
a network of three wide-band, widely separated detec- 
tors was first solved by Giirsel and Tinto [13]. Their 
technique, which is referred to in recent literature as 
the null-stream method, relies on the observation that 
a gravitational-wave burst present in the data of a net- 
work of three wide-band detectors must satisfy a unique 
closure condition. Giirsel and Tinto studied a two- 
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FIG. 1: Geometry of the network and travel times spent by 
a GWB to propagate across a three-detector network (de- 
tectors "H" , "L" , and "V" ) . The locus of constant time delay 
between two detectors form a ring on the sky concentric about 
the baseline between the two sites. For three detectors, these 
rings may intersect in two locations. One is the true source di- 
rection, while the other (S') is its mirror image with respect 
to the geometrical plane passing through the three sites. This 
two-fold ambiguity can be resolved by further considering the 
amplitudes of the responses. For four or more detectors there 
is a unique intersection point of all of the rings. 



parameter family of linear combinations of the three data 
sets, in which the two parameters correspond to the two 
angular coordinates of the hypothesized sky location of 
the source. They showed that when the two parame- 
ters coincide with the true location of the source then 
the gravitational-wave burst is canceled precisely in the 
linear combination. This point is located by applying a 
least-squares minimization (i.e., a test) to the linear 
combination. In [ ] it was also shown that this condition 
holds regardless of the time dependence of the two polar- 
ization waveforms of the burst. This remarkable result 
makes this method very powerful for solving the inverse 
problem since it does not require a priori knowledge of 
the burst waveforms. 

Even with the Giirsel-Tinto technique, lack of knowl- 
edge of the signal waveforms presents a serious im- 
pediment in searches for GWBs. This is because 
gravitational-wave detectors exhibit transient noise fluc- 
tuations, or "glitches" . Without a model for GWB wave- 
forms, it is not obvious how to determine with confidence 
whether a candidate detection is a real gravitational- 
wave burst signal or a "false alarm" due to noise fluc- 
tuations occurring in coincidence in the detectors of the 
network. A confident detection of gravitational-wave 
bursts requires the ability to distinguish them from such 
noise transients. In this context, as pointed out in 
[14, 15, 16], the null-stream method could also be used 



for discriminating gravitational-wave bursts from noise- 
triggered fluctuations affecting the data of the detectors. 
Put simply, uncorrelated noise glitches should not cancel 
in the null stream. In principle, therefore, they can be 
vetoed by setting a threshold on the maximum allowable 
value. In practice, however, a veto test is vulnera- 
ble to effects that prevent precise cancellation of strong 
GWBs, such as calibration errors, and may be ineffective 
for weak glitches, for which there are often sky positions 
which yield x^ ^ ^ V^^ degree of freedom even though 
the glitches are uncorrelated. A threshold which is 
low enough to pass a GWB with poor data calibration 
may also pass a weak glitch. A threshold which is high 
enough to reject the glitch may also reject the GWB. 

In this paper we propose a modified null-stream based 
technique for discriminating GWBs from noise glitches. 
This technique is based on comparing the energy in the 
null stream to that expected if the transients in the de- 
tectors are uncorrelated. This second energy measure, 
which we call the incoherent energy^ provides an effective 
measure of the significance of the x^ test and renders it 
robust against both calibration uncertainties and weak 
glitches. Like the Giirsel-Tinto analysis for determin- 
ing the source direction, our null-stream consistency test 
does not require any a priori knowledge of the GWB (or 
glitch) waveforms. 

Our paper is organized as follows. In Section II we 
show that there exists a very general and elegant proce- 
dure for deriving the null stream for an arbitrary num- 
ber of detectors with colored noise. We then address 
the issue that the null-stream x^ alone can not reliably 
distinguish between a gravitational- wave burst and noise- 
generated glitches. This is done by introducing a com- 
plementary energy measure, the incoherent energy, and 
demonstrating that GWBs and glitches separate in the 
two-dimensional space of null and incoherent energies. 
This allows us to identify (or "veto") noise-generated 
events and hence make the null-stream analysis robust 
against glitches. 

In Section III we discuss the results of the numerical 
simulation of our statistical test applied to the LIGO- 
Virgo 3-detector network. We assume the three inter- 
ferometers to be working at the LIGO design sensitivity, 
and quantify the ability of our method to distinguish true 
GWBs from coincident noise glitches. Although our nu- 
merical implementation is not optimized for the signals 
under consideration, it indicates that gravitational-wave 
bursts observed in each detector with SNRs of about 
10 — 20 can reliably be distinguished from noise glitches 
of similar energy, and that a significant improvement 
over the statistics based on the null-stream x^ alone is 
achieved. Our conclusions and future work plans are pre- 
sented in Section IV. 
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II. ANALYSIS 

Three or more detectors provide redundant measure- 
ments of the two polarization components hx of a 
gravitational wave. It is therefore possible to construct 
linear combinations of the data streams that do not con- 
tain any gravitational- wave component, i.e., that consist 
only of detector noise. In this section we derive these lin- 
ear combinations, known as "null streams," for networks 
containing an arbitrary number of detectors whose noises 
are different and colored. 



A. Conventions 

The conventions used for the notation in this report 
are described in Table I. 

For a plane [ ] gravitational wave incident from a di- 
rection the strain /i+^x at the position fa is related 
to that at some (arbitrary) reference position fo by 



where the time delay Ata{^s) is given by 



Atai^s) -{fo - fa) ■ • 



(1) 



(2) 



We can therefore compare the gravitational- wave signals 
measured by detectors at different locations by shifting 
the time-series data from each detector according to (2), 
provided we know the sky location of the source. 

The time-series signal produced in detector a at fa by 
a gravitational- wave /?^+,x incident from a sky position 

Qs is 



da{t ^ Atai^s)) 



^na{t^Ata{ns))^ (3) 



Here Ua is the stationary background noise of detector 
a as a function of time, calibrated to units of strain, 
and are the antenna response functions [13] for 

detector a for the sky position Qs of the gravitational 
wave source. For brevity, we write h-^^x{t) = /i+,x ^o)- 



B. Null stream construction 

Let us assume for the moment that we know the direc- 
tion Qs to the source. Then we can time-shift the data 
from each detector as in (3), and drop explicit references 
to the time delays Ata and the sky position. Transform- 
ing to the Fourier domain, equation (3) becomes 



(4) 



The Fourier transform and its inverse are defined by 



daif) = 
da{t) = 



+ 00 



dM) e'"'^*d/ 



(5) 



The one-sided strain noise power spectral density Sa{f) 
of the stationary noise Ua is given by 



1 



{na{f)n;{n) = -SapSif - nSaif) . 



(6) 



For present detectors Sa{f) is a strongly varying function 
of frequency. Since it will prove convenient to work with 
white-noise data, without loss of generality we divide the 
strain data at each fr equency by the estimated amplitude 
spectrum ^JSa{J)/^ of the corresponding detector noise. 
The whitened data d^a is then given by 



dwa (/) 



= F+Mf) + F:^ah^ if) + hUf) , (7) 



where the n^ait) are unit Gaussian noise processes and 
F^^^ are the noise-weighted antenna responses 



_ F+'-{hs) 



(8) 



The F^^ contain all of the information on the detector 
sensitivity, both as functions of frequency and source sky 
position. 

For a network of D detectors, equation (7) can be writ- 
ten in the equivalent matrix form 



dwi 




dw2 




_ dwD _ 





^w2 ^w2 

















nw2 




hx 








_ ^wD _ 



(9) 



or 



(10) 



where we use boldface to denote vectors and matrices. 
Here 



h 



(11) 



and the matrix is defined as 



^wl ^wl 
^w2 ^w2 



77+ 77X 



(12) 
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D 




Number of detectors in the network. 


a,(3 G [1,.. 




Index specifying detector. 


X 




(Boldface) A vector or matrix on the space of detectors. 






Matrix transpose of X. 


N 




Number of data samples from each detector being analyzed. 


j,/cG[0,... 




1] Index specifying time or frequency sample. 






"Plus" and "cross" polarization waveforms of the gravitational wave. 






Sky position of the gravitational- wave source. 






Trial sky position. 


d, n, F 




Data, noise, and antenna responses in strain. 


dyM 5 5 -F'w 




Noise- weighted (whitened) data, noise, and antenna responses. 



TABLE I: Notation conventions for commonly used quantities in this paper. 



This form makes it clear that, regardless of the functional 
form of /i+,x, the gravitational- wave burst can only con- 
tribute to the network output along the directions 
and . The construction of null streams is thus obvi- 
ous: we simply project the data d orthogonally to these 
directions. Formally, we select a new orthonormal Carte- 
sian coordinate basis for the space of d in which vec- 
tors cd-i and bd span F^ and F^ . The remaining basis 
vectors ei, ... ^€0-2 are then orthogonal to and F^ , 



F+ 







e{l,...,D-2}. (13) 



We say that the ei=i^...^i:)_2 form an orthonormal basis for 
the null space of F^^ [33], hence the term "null stream 
formalism." We then construct di {D — 2) x D matrix 
A whose rows are the components of this orthonormal 
basis, 



-D-2 



(14) 



By construction A is orthogonal to and F^ , so 

AF^ = 0. (15) 

We obtain the null streams z by applying A to the net- 
work data vector: 



AF^h- 
An^ . 



An^ 



(16) 



The two independent strain components h-^^x are can- 
celed out as a consequence of the definition (15) of the 
null space, making each Za a null stream. 

For example, the three-detector case is particularly 
simple. In this case the matrix A is equal to 



A = 



F+ X F> 

w w 

\F+ X Fi 



(17) 



This is illustrated schematically in Figure 2. In the gen- 
eral case A can be obtained either via singular value de- 
composition [17] or by explicitly constructing the associ- 
ated projection operator (see Appendix A for details). 

Note that A is a function of the sky position (through 
F+'X (!],)) and frequency (through Sa{f))^ The GWB 
will only be canceled when A is evaluated for the correct 
source location, since generally 



(18) 



In the next subsection we will discuss how to deal with 
the case where the source location is not known. 

In the preceding discussion we have assumed implicitly 
that F^ and F^ are independent. In the general case 
the number of independent null streams is D — where 
r is the number of independent columns of F, i.e. 



rank(Fw) . 



(19) 



There are two cases 



1. If at least one detector in the network has a dif- 
ferent alignment from the others then F^ and F^ 
are independent and rank(Fw) =2. In this case 
there are D — 2 null streams, and the method is 
applicable to networks of D > 3 detectors. 

2. If all detectors in the network are aligned then 
F^ (X F^ and rank(i^w) = 1. In this case there are 
D — 1 null streams, and the method is applicable 
to networks of D > 2 detectors. 

Although in this paper we will concentrate on case 1, i.e. 
three or more non-aligned interferometers, we emphasize 
that our method works also with aligned detectors. 



C. Null-stream analysis 

Since the are orthonormal by construction, it follows 
that 



AA' =1 



(D-2)x(D-2) • 



(20) 
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FIG. 2: Geometry of the null stream construction for the 3- 
detector case. The null stream is obtained by projecting the 
data along the vector A, which is orthogonal to and . 
For D non-aligned detectors A has D — 2 dimensions. 

This implies that each is a Gaussian random process of 
unity variance and is uncorrelated with for dl\ j3 ^ a. 
This is the main advantage of the null-stream formalism 
over other techniques: the noise distribution of the pro- 
jected network data is known a priori, regardless of the 
form ofhj^^x, under the assumption that a GWB from a 
particular direction is present. This allows us to perform 
statistically significant tests of the network data, in par- 
ticular to test the hypothesis that a GWB from a given 
direction is present. 

Since data from gravitational-wave detectors are sam- 
pled and digitized, in what follows we will consistently 
use discrete notation in our analysis of the statistics of 
the null stream. Our conventions for discretely sampled 
data are as follows: the Fourier-transform pair becomes 

N-l 

x[k\ = ^x[j]e-^2^^'^/^, 

1 ^-^ 

^[^'] = ]v E^W^^'"''^"^' (21) 

k=0 

where is the number of data points in the time domain. 
Denoting the sampling rate by /g, we can convert from 
continuous to discrete notation using x{t) x{f) 

f-'x[k]. Jdt^ f-' E,, / # - fsN-' Efc, ^{t - n - 



fsSjj', and S{f - f) Nf^ ^hk'- For example, the 
one-sided strain noise power spectrum Sa[k] is 

N 

{K[k]Mk']) = Y^-f^^kk'Saik] . (22) 

We will whiten the data by applying a zero-phase whiten- 
ing filter [18, 19], and our normalization convention for 
whitened data is 

(^wa[^]^w/3[^1) = Saphk' • (23) 

The total energy in the null streams is 

D-r N-l 

^nun^EEl^'"Wl'- (24) 

a=l k=0 

Using (20) and (23) it follows that at the true source 
position 2£^nun is x^-distributed with 2N{D — r) degrees 
of freedom. In this case the expectation value of the null 
energy and its variance are both N{D — r). 

Although our considerations so far have assumed the 
sky position of the GWB source to be known a priori^ 
in practice this may not be the case. Since we know that 
at the correct source location the null energy will not con- 
tain any contribution from the signal, a straightforward 
procedure is to scan over a grid of sky positions in search 
of the minimum of the null energy. Giirsel and Tinto [13] 
used time-delay estimates to limit their search to two 
possible regions of the sky (the regions around points S 
and S' in Figure 1). This approach may fail when the 
duration of the signals (or the timing uncertainty) is of 
the same order as the light travel time between detectors, 
thus necessitating an all-sky search. Since the numerical 
analysis in [ ] implies that the characteristic angular 
width of the minimum of the null energy in a neighbor- 
hood of the source location for a GWB with a central fre- 
quency of about 100 Hz is equal to approximately 10~^ 
steradian for high SNRs, it follows that an all-sky search 
should be performed over a sky grid containing more than 
10^ resolvable directions. (Our numerical tests use a grid 
containing 10^ points.) In either case, for each trial direc- 
tion one postulates the presence of a gravitational-wave 
signal, forms a linear combination of the detectors that 
is orthogonal to that postulated direction, and x^-tests 
this null stream for excess energy. If there exists a par- 
ticular direction for which there is no excess energy in 
the null stream, the data is regarded as consistent with 
the hypothesis that a gravitational- wave burst is present 
and incoming from the inferred direction. If, on the other 
hand, £^nuii is inconsistent with a x^ distribution, then 
one rejects the hypothesis that a GWB is present incom- 
ing from that direction. The best estimate of the source 
direction is taken as the direction with minimum x^- 

Although the null-stream method does not require 
knowledge of the two GWB waveforms for its imple- 
mentation, once the source location Qs has been iden- 
tified it is straightforward to reconstruct h-^^x from the 
data themselves [ ] (if the detectors are all aligned then 
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only one of the polarizations can be reconstructed.) The 
minimum-variance estimate of the two waveforms for a 
network containing an arbitrary number of non-aligned 
detectors is given in Appendix B. 

The null-stream combination of the data from a three- 
detector network and the resulting test were first de- 
rived (in a different way) by Giirsel and Tinto [ ] for 
detectors whose noises are white. They also implemented 
a near-optimal filtering procedure to account for colored 
noise. Our approach builds on this by generalizing and 
simplifying the derivation of the null stream and wave- 
form reconstruction to networks containing an arbitrary 
number of detectors with different colored noises. One 
can also show that the null-stream procedure is formally 
equivalent to the maximum-likelihood analyses presented 
in [6, 20], though we leave the demonstration to a future 
paper. 



D. Distinguishing GWBs from Noise Transients 

The numerical analysis performed by Giirsel and Tinto 
in [13] was not aimed at checking whether the null energy 
estimator could distinguish GWBs from noise-induced 
glitches. In fact, an analysis based purely on the null 
stream runs into difficulty when applied to data contain- 
ing noise transients. Strong uncorrelated glitches gener- 
ally will not cancel in the null stream combination be- 
cause they are not correlated in amplitude and phase 
in a way consistent with a GWB, implying a > 1 
per degree of freedom. However, a null-stream analysis 
of a real GWB may also produce x^ > 1 per degree of 
freedom, due to imperfect cancellation of the GWB in 
the null stream. This may happens for various reasons, 
such as the use of a discrete sky grid, inaccurate calibra- 
tion of the data, or imperfect whitening of non-stationary 
data. Thus, an analysis based purely on the null stream 
would be forced to either reject both glitches and GWBs 
or accept both. Further, counter to one's intuition, this 
problem could get worse with stronger signals. 

Glitches can also fool a null-stream analysis in 3- 
detector networks when the transient is weak in at least 
one detector. This is because the non-observation of a 
signal by one detector a is always consistent with a GWB 
incident from a direction and polarization to which the 
detector responses F^^ and F^^ are sufficiently small. 
For sky positions with F^^ ^ F^^ ~ the null stream 
projection matrix A reduces to 



1/37 



wa 1 wa 



0. 



(25) 



That is, the null stream for this sky position reduces to 
the detector in which there is no transient. This gives 
a x^ per degree of freedom of order unity regardless of 
whether the transient in the other two detectors is due 
to a GWB or noise. (Equivalently, when only two detec- 
tors observe a signal, it is always possible to find h-^^x 
that fit the output of these two detectors.) Thus double- 
coincident glitches will always pass a x^ test for certain 



areas on the sky. And while networks containing four or 
more detectors will be less affected by this problem [21] 
because the size of the region of the sky producing two or 
more simultaneous, below-threshold responses is smaller 
than that for a single response, it will not be null. 

In what follows we propose a simple way to make null- 
stream analyses robust against glitches by comparing the 
amount of energy in the null streams to that expected 
if the transients are uncorrelated. Let us consider how 
the null stream energy £^nuii depends on the individual 
detector data streams (i^a- By defining the matrix 



Q = A'A, 

we may write (24) in the convenient form 

N-l D D 
/c=0 a^l /3=1 



(26) 



(27) 



N-l 



k=0 



Qu Ql2 • • • QlD 

Q21 Q22 Q2D 

QdI Qd2 • • • QdD 



dw2 



vD 



(28) 



Note that the null energy contains contributions from 
both cross-correlation {d^^^d^p) and auto-correlation 

{d^^d^^) terms. If the signals in the various detectors 
are independent (as one might expect for noise glitches), 
then the expectation value of the cross-correlation terms 
in (28) will be small compared to that of the auto- 
correlation terms. In this case the expectation value of 
FnnW is just the sum of the diagonal terms in (28): 



N-l D 



mean(£;nuii(^)) — > ^ ^ Qaa[k,'^]\d^a[ 



(29) 



Here the mean is an ensemble average over noise instan- 
tiations. This observation motivates the use of a new 
energy measure, the incoherent energy E'inc, defined as 
the autocorrelation contribution to the null energy: 



N-l D 



E,nc{^) = ^ ^Q«c.[^,^]|^iwa[^]|'. (30) 



If the transient signals in the various detectors are not 
correlated, then we expect the following approximate 
equality to hold: 



mean(£;nuii) - mean(£;inc) • 



(31) 



If instead the signals in the detectors are correlated, as 
will be the case when a GWB is present in the data, 
then at the correct sky location the GWB contributions 
cancel in the null energy, and the following inequality 
should hold: 



mean(£^nuii) < mean(£^inc) • 



(32) 
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Thus, a distinguishing feature of a GWB is that a sig- 
nificant fraction of the energy in the individual detector 
data streams is canceled in the null stream. 

In our simulations we test two simple measures of the 
degree to which £^nuii and E[nc show the behavior ex- 
pected of a GWB. These are the quantities £^nuii — ^inc 
and (^^nuii - Einc)/Einc = £^nuii/^inc " 1, which represent 
respectively the amount of "correlated energy" and the 
ratio of "correlated energy" to "uncorrelated energy" in 
the data. 

An example is shown in Figure 3. This figure shows 
the incoherent energy versus null energy for a GWB and 
a glitch of the same amplitude, evaluated for approxi- 
mately 10^ uniformly distributed sky locations (see III for 
details). We note that for both the GWB and the glitch 
there are sky positions for which the null energy is ap- 
proximately unity, so that projection onto the null energy 
alone does not distinguish this glitch and GWB. However, 
glitches and GWBs scatter differently in terms of the two 
energy measures E'nuii and £^inc- This suggests a modified 
procedure for distinguishing GWBs from glitches: scan 
over the sky and look for directions for which £^nun is 
significantly smaller than £^inc- If there exists a sky di- 
rection for which £^nuii is sufficiently small compared to 
^-inc, we conclude that the transient could indeed be a 
GWB. If instead there is no direction for which £^nuii is 
sufficiently small compared to Ei^c , we conclude that the 
transient is not a GWB. 

Strong GWBs that are not precisely canceled in the 
null stream (due to calibration errors, for example) and 
have £^nun > N{D — r) will still pass this test because 
^nuii < Einc Double-coincident glitches will fail even 
though £^nuii — N{D — r) because they have £^nun — 
^inc- Put another way, the incoherent energy provides 
a natural cutoff in the significance of a P^i" degree 
of freedom measurement. Intuitively, if the null stream 
cancels out some large fraction of the excess incoherent 
energy, then we expect the event to be a gravitational 
wave, even if > 1 per degree of freedom. The failure to 
cancel a significant portion of the incoherent energy will 
eliminate glitches even if — 1 per degree of freedom. 

To get more insight into the behavior and usefulness of 
these energy measures, let us consider in more detail their 
expectation values over noise instantiations. Allowing for 
the case in which the transient is not a GWB, we write 



[k] = n^a[k] + gwa[k] 



(33) 



where g^a denotes the noise-weighted transient as seen 
in detector a. For example, if the transient is a GWB 
from the direction (^g, and the tested sky position is 1], 
then 



9wa [k] 



X e 



227r/s/e/Ar(At«(^)-At«(^s)) 



(34) 



incorrect sky location Q instead of the true but unknown 
sky location Qg. 

For simplicity, let us restrict ourselves to the case of 
detectors with equal noise spectra, for which Q is inde- 
pendent of frequency. Then the noise-weighted transient 
signal g^a appears in the null and incoherent energies in 
the combination 



N-l 



AT-l 

S 'n 



,(35) 



which is the noise-weighted cross-correlation of the sig- 
nals in detectors a and p. The diagonal terms are 



N-l 



N-l 



2 

Paa 



\Uk]? 



(36) 



This autocorrelation term is simply the squared signal-to- 
noise ratio of an optimal matched filter for the transient 
in detector a, as indicated by the second line of (36). 

As an example, let us consider the special case of a 
linearly polarized GWB (e.g., with hx =0), with trial 
sky position Q and true source position Qg- In this case 
the cross-SNR is equal to 



2 

TV 



N-l 



^Sa[k]S0[k] 



COS'H al3{^,^s) (37) 



(only the real part of p'^^ contributes to the energies), 
where the phase error is 



27r^(At„(f}) 



(38) 



The phase term accounts for time shifting based on the 



The cross-SNR for a 7^ ^ is typically of the same 
order of magnitude as p^^,, but it is positive or negative 
depending on the timing error between pairs of detec- 
tors. As we test different positions on the sky, the tim- 
ing errors change and the GWB contributions from the 
different detectors move in and out of phase. This will 
produce interference fringes in maps of E^uW^ £^nuii — ^inc, 
and £^nuii/^inc- The location and spacing of these fringes 
are determined by the dominant frequency of the signal, 
the true sky position of the source, and the detector ge- 
ometries and relative locations. 

In terms of the SNRs (35), (36), the expectation values 
of the lowest moments of E'nuii, ^inc, and £^nuii — Ei^c are 

D D 

mean(Enuii) = iV(I^ - r) + ^ ^ Q,;3p'/3 , (39) 

a=l (3=1 
D D 

var(Enuii) = A^(I^-r) + 2^^Q,;3P^^,(40) 

a=l /3=1 
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FIG. 3: Scatter plot of the null and incoherent energies for a simulated GWB (left) and a simulated glitch (right) seen in a 
network consisting of the LIGO-Hanford and LIGO-Livingston 4 km detectors and Virgo. In both cases the signal was scaled 
to rms SNR of 20 over the three detectors. Each point represents one trial sky position; approximately lO'^ sky positions in 
a uniform grid were tested. The waveforms and detector network used are discussed in detail in Section III. The GWB and 
glitch signals have the same signal-to-noise ratios and time delays in the individual detectors, and so are indistinguishable to 
incoherent tests. Note that for both the GWB and the glitch there are sky positions for which the null energy is consistent 
with noise (^nuii — N{D — r) — 60 for these simulations). However, for the GWB there are also sky positions with £^nuii ^ ^inc 
(points above the diagonal); these are due to the fact that the GWB signal is correlated between the detectors. The glitch 
signal does not access this portion of the (^nuii, ^inc) space. This observation is the basis of our consistency test; we scan over 
the sky and look for directions where ^nuii < ^inc- The true source location is indicated by the circle on the GWB plot. Also 
shown on the GWB plot are the tested sky positions which have the smallest values of ^nuii — ^inc, ^nuii/^inc, and ^nuii- The 
sky maps for these same simulations are shown in Figure 4. 



D 

mean(£;inc) = N{D - r) ^ ^Q^apla ^ (41) 

D D 

var(i?i„e) = iv^gL + 2^gLpL, (42) 

D D 

mean(^null - ^inc) = ^^{'^ - ^a/3)Qa/3pl(3 ^ (43) 
a=l/3=l 
D D 

Var(^null - ^inc) = A^^^(l-^c./3)Q'/3 + 

a=l /3=1 

D D D 

+ 2 E E E(l - ^«/5)(l - ^a^)QapQa^pl^ • (44) 
a=l /3=1 7=1 

We note that the signal enters E^nc only through its SNR 
p^Q, in the individual detectors; E^nc does not depend on 
the structure of the transient signal. As a result, varia- 
tions in E^inc reflect variations in the network sensitivity 
due to noise-weighted geometrical factors (the Q) and do 
not contain significant information on the signal. By con- 
trast, £^nun — Einc contains only cross terms, and shows 
the interference of the signals measured by the different 
detectors. 

For example, Figure 4 shows sky maps of the £^nuib 

^inc, and £^null/^inc [notC that £^null/^inc = (£^null " 



Einc)/ Eiric + 1] for the same GWB and glitch signals used 
in Figure 3. The null energy maps for the GWB and the 
glitch are very similar. The GWB and glitch are con- 
structed to have the same relative time delays and SNRs 
in the various detectors. As a result, the incoherent en- 
ergy maps are virtually identical for the GWB and the 
glitch. Removing this signal-independent structure from 
the null energy makes the signal-dependent structure in 
the sky maps clearer. In particular, the plot of ^nuii/^inc 
for the GWB shows sharp interference fringes orthogonal 
to the Hanford-Livingston and Livingston- Virgo base- 
lines (the signal was strongest in Livingston and Virgo 
in this simulation). The sky location of the GWB signal 
lies on one of the two intersection points of these inter- 
ference fringes, so they can be used to locate the source. 
Such sharp features are not present in the corresponding 
sky map for the glitch, since the glitch waveforms are not 
strongly correlated. 

We note from (39)- (44) that fractional fluctuations 
in the energies scale as A^~^/^ for weak signals. Since 
the signal-dependent structure scales as p^^, the basic 
limit of purely local measurements of sky maps scales as 
p\^lN^I^ . This is the same scaling as in excess-power 
searches [22]. 

Also, we note that if we restrict our analysis to a 
frequency range / = /c/^/TV G [/min, /max], then TV in 
(39)- (44) becomes the number of frequency bins actually 
summed over, and p^ (35), (36) is to be computed over 
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FIG. 4: Sample sky maps, normalized to E/(N{D — r)), for the same GWB and glitch signals used in Figure 3. The plots 
on the left are for a GWB with rms SNR of 20 in the three detectors. Those on the right are for a glitch with the same 
relative time delays and signal energies in each detector as the GWB. (Note that these two events are indistinguishable to 
an incoherent analysis.) The network consists of the LIGO-Hanford (H) and LIGO-Livingston (L) 4 km detectors and Virgo 
(V). The null energy map (top) shows interference fringes due to the transient signal, as well as structure due to the network 
geometry. They are very similar for the GWB and the glitch. The incoherent energy map (middle) is constructed from the 
auto-correlations of the individual detector data streams and reflects variations in the network sensitivity over the sky. It is 
virtually identical for the GWB and the glitch because the two events have the same relative time delays and SNRs. Dividing 
out the incoherent energy from the null energy (bottom) removes this structure associated with the network geometry, making 
the signal-dependent structure clearer. The GWB map shows sharp interference fringes (blue and red rings) where the time 
delays along the H-L, L-V, and H-V baselines match those of the source location. There is little of such structure in the glitch 
map. 



10 



the same (positive) frequency range. 

Finally, we point out that in the high SNR limit the 
minimum of the null energy £^nuii occurs at the source 
location [13], which is useful for solving the inverse 
problem for bursts. The incoherent energy Ei^c is not 
an extremum at the source location, so measures like 
£^nuii/^inc cLud E^nuii ~ ^inc cauuot dircctly resolve the 
source location. As noted above, however, the clearer 
signal-dependent structure in £^nuii — ^inc and £^nuii/^inc 
compared to £^nuii may be useful for this purpose. In 
any case, our consistency test can improve the detection 
confidence of a gravitational- wave burst, and should be 
regarded as an essential first step in solving the inverse 
problem. 



E. Sensitivity of the method to data calibration 

The derivation of the null-stream combinations de- 
scribed in subsection JIB assumed the data from the 
interferometers to be perfectly calibrated. This means 
that properly modeled transfer functions of the interfer- 
ometer responses (which depend on various parameters 
associated with the characteristics of the cavities of the 
interferometers) are applied to the raw data in order to 
obtain strain measurements. In practice, however, the 
parameters describing these transfer functions are known 
with finite accuracy, in turn preventing the null-stream 
combinations from exactly canceling the GWB signal at 
the correct source location. In order to quantify the ef- 
fect on our method we should note that the magnitude of 
the residual signal in the null-stream combinations will 
be proportional to the accuracy by which the calibration 
parameters are known. For the LIGO detectors the cali- 
bration parameters are known to better than ten percent 
[23]. Preliminary studies indicate our consistency test is 
robust against calibration errors of this size; these effects 
will be studied in more detail in a future article. 



F. Implementation 

Our null-stream based analysis has been implemented 
as a publicly available Matlab package, "Xpipeline" 
[24]. For a specific detector network and event time, 
Xpipeline reads the appropriate data from frame 
files [25] (the standard format for storing data from 
gravitational- wave detectors), optionally injects GWB 
signals and/or glitches, whitens the data, computes the 
null stream coefficient matrix A for each specified sky di- 
rection and frequency, computes the time shifts for each 
direction, steps through data in overlapping blocks of 
user-specified duration, time shifts the data to the near- 
est sample, Fourier transforms it, completes the time 
shift with a phase rotation, forms the null stream in 
the frequency domain, sums the power in user-specified 
frequency bands, and records the null and incoherent 
energies for each time-frequency band and direction. 



Xpipeline runs in approximately 1/lOOth real time. For 
example, the analysis of the 10^ simulated events used to 
produce Figures 6-11 took approximately 16 hours on 4 
Intel Pentium 4, 2.66 GHz computers. This makes our 
null-stream-based consistency test feasible as a follow-up 
test in GWB searches. 



III. SIMULATIONS 
A. Network and signal types 

To test the efficacy of our statistical test in discriminat- 
ing GWBs from noise glitches in the (^^nuih^inc) space, 
we need to select a detector network, a population of 
GWBs, and a population of glitches. 

We elect to simulate a network consisting of the Han- 
ford and Livingston 4 km interferometers ( "H" and "L" ) 
and a third identical instrument at the Virgo site ( "V" ) . 
For the sake of simplicity we neglect both the Han- 
ford 2 km interferometer and the differences between the 
LIGO and Virgo design sensitivities. The locations fa 
and orientations (which determine F+, F^) of the in- 
terferometers are taken from [22, 26]. The calibrated 
stationary background noise ria [j] for H and L are taken 
from a standard 24 hour reference simulation [27]. The 
background noise for V is taken from the Hanford sim- 
ulation with a 2 second time shift, which is much larger 
than the time scales of the signals used in this analysis. 

The next step is to select the GWB and glitch wave- 
forms. To simulate a glitch, we need three waveforms 
(one for each detector) that are not strongly correlated. 
To simulate GWBs we wish to use waveforms that are 
motivated by astrophysical considerations. Our consis- 
tency test is ultimately based on the fact that for a GWB 
the signal seen in each detector is correlated in a particu- 
lar way, whereas for a glitch the signals generally will not 
be correlated. In order to show that our test does not rely 
on any fundamental difference between "typical" GWB 
waveforms and "typical" glitch waveforms, we elect to 
use the same set of waveforms for simulating GWBs and 
glitches. We select three representative waveforms from 
the Dimmelmeier- Font-Mueller ("DFM") catalog [9] of 
Type II core-collapse supernovae. Specifically, we choose 
the A1B3G3 "regular collapse" waveform, the A1B3G5 
"rapid collapse" waveform, and the A3B4G2 "multiple 
bounce" waveform. Figure 5 shows the time-series and 
power spectra of these waveforms. As we shall see below, 
these three waveforms have moderately low but nonzero 
cross-correlations. 

To simulate a gravitational-wave signal, one of the 
three waveforms is randomly selected and added into the 
data stream from each interferometer, with time delays 
and scaling appropriate for some choice of polarization 
angle and location on the sky. To simulate a glitch we fol- 
low the same procedure, except that a different waveform 
is selected for each detector. The scaling and time delays 
proceed as for a GWB. This population of glitches has 
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FIG. 5: Dimmelmeier - Font - Mueller ("DFM") waveforms used to simulate GWBs and glitches in our analysis. Left plot: 
time series waveforms scaled to SNR = 100 for optimal orientation. (The linear trend present in the original waveforms has 
been removed.) Right plot: spectra with a 1 Hz resolution; shown for comparison is the LIGO design noise curve used for both 
LIGO and Virgo detectors. 



the property that it would pass any incoherent test, as 
the arrival times, power distribution, and even frequency 
bands are consistent with those of true gravitational- wave 
signals. We stress that, though not a realistic glitch pop- 
ulation, this provides us with examples of the kind of 
pathological glitches that cannot be dismissed by per- 
detector methods. 



sky with a maximum angular mismatch defined by the 
longest detector baseline and the maximum frequency. 
As with all template placement problems we have some 
freedom in how to produce this set, and chose a some- 
what sub-optimal but simple set of directions - a grid of 
approximately 10^ points uniformly distributed in 9 and 
0sin^. 



B. Analysis parameters 



C. Waveform normalization 



The main tunable parameters in our analysis are the 
time and frequency bands and the sky positions to test. 
Ideally the integration time and frequency band should 
be matched to the signal being tested, to maximize the 
signal-to-noise ratios p^^, and minimize the background 
noise contribution N{D — r). We choose an integration 
length of 1/16 s, which is the smallest power of 2 larger 
than the durations of the three sample waveforms tested. 
We overlap consecutive data segments by 50% to min- 
imize the loss of signal-to-noise when a signal overlaps 
the edge of a data segment. We use a single frequency 
band [34] of [64, 1024] Hz. The upper frequency limit is 
set by the highest frequency at which our target signals 
have significant power compared to the noise curve; see 
Figure 5. The lower limit is set at 64 Hz because the 
actual noise level in current detectors is larger than the 
design noise below this cutoff [28]. 

Once the locations of the detectors and the upper 
bound on frequencies involved in the analysis are known, 
a set of directions {ft} covering the sky may be produced. 
Both the projection matrix Q (26) and the time delays 
Ata (2) vary with angle, but the effect of At^ on the 
cross-correlation terms occurs on a smaller angular scale 
for most of the sky. A simple criteria then is to cover the 



The DFM waveforms are linearly polarized, which 
means that the two polarizations are linearly dependent 
and so can be written in the form 



(45) 
(46) 



where ips is the polarization angle. As a result, the strain 
signal ga in detector a is 



a{t)= (^F^{hs)co82^s^F^{hs)sin2^s)h{t). (47) 



9> 



These waveforms are de-trended and normalized against 
the interferometer noise curve so that p^^, = 1 (36) for 
optimal orientation (the case cos 2ips + sin 2ips = 
1, so that ga{t) = h{t)); i.e., we define the normalization 
of h{t) so that 



AT-l ,7 



mi 

.^0 TS{k] 



E 



1. 



(48) 



(Recall that our detectors have identical noise spectra, so 
Sa[k] = Sf3[k] = S[k].) The cross-correlations (35) of the 
waveforms depend on their relative time or phase shift; 
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with this normahzation the maximum cross-correlations 
over ah shifts for co-ahgned detectors are 



max^ {pip} = 




(A1B3G3 
(A1B3G3 
(A1B3G5 



A1B3G5) 
A3B4G2) 
A3B4G2) 



(49) 



For comparison, typical cross-correlation values for Gaus- 
sian noise in our time- frequency band are 0.15 — 0.2. 

To simulate a gravitational wave signal, one of the 
three waveforms A1B3G3, A1B3G5, or A3B4G2 is ran- 
domly selected. Uniformly distributed random direction 
fls and polarization tps angles are chosen. For each detec- 
tor a, the discrete catalog waveform h{t) is time-shifted 
by Atai^s) (2) and resampled to match Ualj]- The wave- 
form is then scaled by the antenna response as in (47) to 
give ga- 

To characterize the efficacy of our consistency test as 
a function of the signal strength, we choose to simulate 
populations of candidates with the same measured signal- 
to-noise ratios; i.e., the signals are scaled so as to deliver 
a fixed total SNR to the network. This simulates candi- 
dates near a detection threshold, rather than (for exam- 
ple) a physical population of standard candles. With the 
normalization (48) and (47) one finds 



D 



D N-l 



D 



= {Fa{^s)cos2^s^F^{hs)siri2^sy .(50) 



To fix the SNR in the detectors to some rms value Prms, 
we apply the further normalization 



9a{t) 



D 

\ i^(^/3^(^.)cos2V;. + F^^(l^.)sin 



With this scaling we find 



D 



D ^ 



(51) 



(52) 



To simulate a population of glitches, the same process 
was followed with the sole exception that a different DFM 
waveform was selected for each detector. We applied the 
same scaling and time delays as for a GWB. This popu- 
lation of glitches has the property that it would pass any 
incoherent test, as its arrival times and power distribu- 
tion are consistent with those of true gravitational- wave 
signals. We reiterate that, though not a realistic glitch 
population, this provides us with examples of the kind 
of pathological glitches that cannot be dismissed by per- 
detector methods. 



D. Analysis Procedure 

After the signal has been added to the background 
noise, the data are whitened to produce (iwa[j]- The 
whitening algorithm is trained on a 16 second block of 
data that does not include the signal. With a known set 
of trial sky positions {(7} and measured power spectra 
Sa{f)^ Q can be computed for each direction and resolv- 
able frequency. 

For each direction on the sky, overlapping segments of 
data are considered sequentially, the length of the seg- 
ments depending on the time scales of the signal under 
consideration (here chosen as 1/16 s). To perform the 
time shifts, the segments are extracted from the clos- 
est integer samples in the time domain, and then trans- 
formed to the Fourier domain where the remaining part 
of the time shift is performed by applying phase shifts. 
The null and incoherent energies are then computed and 
recorded for that direction on the sky. 

For each simulation we select the two directions for 
which the transient shows the most correlation according 
to two different criteria, as well as the direction with the 
minimum null energy: 

1. min(£^nuii — £^inc)* Measures the linear distance 
away from the diagonal in a scatter plot of ^inc 
vs. £^nuii (see Figure 3). Physically it represents 
the largest amount of energy canceled in forming 
the null stream. 

2. min(£^nuii/£^inc)- Measures the angular distance 
away from the diagonal in a scatter plot of Ei^c 
vs. ^nuii (see Figure 3). Physically it represents 
the largest fraction of energy canceled in forming 
the null stream. 

3. min(£^nun)* Calculated for comparison against the 
other two statistics, and for estimating the sky lo- 
cation. 

Figures 6, 7, 8 show scatter plots of Ei^c vs. £^nuii for 
the sky locations picked by these three criteria for 10^ 
simulated signals. The signal population consists of 10^ 
GWBs and 10^ glitches at each of 5 different signal-to- 
noise ratios, prms = 5, 10, 20, 50, 100. 

Note that there is a significant difference in the distri- 
butions of signals and glitches using any of the measures 
min(£;nuii - £^inc), min(£;nuii/£^inc), or min(£;nuii). For 
both GWBs and glitches there are directions on the sky 
for which the null energy is low. However, for glitches 
these occur exclusively at low incoherent energies. For 
signals there exist directions with low null energy at 
larger incoherent energies. This restriction is strongest 
when measuring correlations using min(£^nuii/^inc) (Fig- 
ure 6). For min(£^nuii — ^inc) (Figure 7) there are a small 
but noticeable fraction of GWBs that overlap the glitch 
population and so are indistinguishable from glitches. By 
using instead min(£^nuii) only (Figure 8) a significant frac- 
tion of the glitches have a sky location with low null en- 
ergy, but only at low incoherent energy as well. 
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FIG. 6: Scatter plot of the null and incoherent energies for 
the most correlated direction on the sky, defined as the direc- 
tion of min(£^nuii/^inc) (i.e., the upper left limit of Fig. 3), 
for glitches and GWBs of various signal-to-noise ratios. Note 
that as the signal-to-noise ratio increases (higher Einc) the 
GWB and glitch distributions separate, with the GWBs re- 
maining at relatively low null energies and the glitches hav- 
ing comparable null and incoherent energies. The clumping 
is due to the 5 distinct rms SNRs used in our simulations: 
5, 10, 20, 50, 100. The null energies for some of the GWB sim- 
ulations are greater than that expected (60) from the number 
of degrees of freedom because we are selecting sky positions 
from the minimum of £^nuii/^inc, not from the minimum of 

^null- 



FIG. 7: Scatter plot of the null and incoherent energies for 
the most correlated direction on the sky, defined as the direc- 
tion of max(£^inc — ^nuii) (i-c, the upper left limit of Fig. 3), 
for glitches and GWBs of various signal-to-noise ratios. Note 
that as the signal-to- noise ratio increases (higher Einc) the 
GWB and glitch distributions separate, with the GWBs re- 
maining at relatively low null energies and the glitches having 
comparable null and incoherent energies. The clumping into 
horizontal bands is due to the 5 distinct rms SNRs used in 
our simulations: 5, 10, 20, 50, 100. The null energies for many 
of the GWB simulations are greater than that expected (60) 
from the number of degrees of freedom because we are select- 
ing sky positions from the minimum of ^nuii — ^inc, not from 
the minimum of ^nuii- 



As the SNR increases the populations become dis- 
tinct. We can see this in Figures 9-11, which show the 
Receiver- Operator Characteristic (ROC) curves for the 
performance of our various statistics in discriminating 
the two populations. At total energies corresponding to 
10 — 20 or greater we can detect most of a pop- 
ulation of gravitational waves and reject essentially all 
of a population of semi-correlated glitches. The rejected 
gravitational waves are those that are weak in at least 
one detector, or equivalently produce little correlation. 

For prms = 5 or lower the glitch and GWB populations 
are not distinct. This is because the null stream enforces 
waveform consistency only when there is excess energy 
to suppress. 



IV. CONCLUSIONS AND FUTURE 
DIRECTIONS 

In this article we have introduced an extension of the 
Giirsel-Tinto null-stream technique that allows one to 



make robust tests to distinguish between GWB signals 
and coincident noise glitches. This technique is based 
on comparing the energy in the null stream to that ex- 
pected if the signals in the various detectors are uncorre- 
lated, and does not require any a priori knowledge of the 
GWB or glitch waveforms. We applied this technique to 
the case of the LIGO- Virgo 3-detector network at design 
sensitivity, and quantified the ability of three different 
measures of correlation to distinguish true GWBs from 
coincident noise glitches. For the best-performing mea- 
sure, the ratio of null energy to incoherent energy, we 
found that gravitational-wave bursts of SNR 10 — 20 or 
greater can be distinguished from glitches of comparable 
SNRs that are injected in the data with the same time de- 
lays but are different in the three detectors' data. For ex- 
ample, with the GWB and glitch populations tested, we 
found that 90% of glitches can be rejected while accepting 
94% of SNR 20 GWBs and 76% of SNR 10 GWBs. Fur- 
thermore, we stress that the glitch population tested was 
pathological in the sense that they were constructed to 
have time delays and amplitudes consistent with a GWB. 
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FIG. 8: Scatter plot of the null and incoherent energies for 
the best-fit direction on the sky, defined as the direction of 
min(^nuii) (i-e., the left-most point in Figure 3), for glitches 
and GWBs of various signal-to-noise ratios. Note that as 
the signal- to- noise ratio increases (higher ^inc) the GWB and 
glitch distributions separate, with the GWBs remaining at 
relatively low null energies and the glitches having comparable 
null and incoherent energies. The clumping is due to the 5 
distinct rms SNRs used in our simulations: 5, 10, 20, 50, 100. 
Note that a large fraction of the glitch signals produce low 
null energies for some sky positions. This means that GWBs 
and glitches are not distinguishable using only the null energy. 



Hence, the performance of the consistency test may be 
even better with real detector data. This consistency test 
is therefore a promising technique for rejecting noise co- 
incidences and increasing detection confidence in GWB 
searches. 

The development of coherent analysis techniques for 
GWB detection is still at an early stage, and much fur- 
ther research can be done. In this section we briefly out- 
line some of the directions of current and future work. 
These can be divided roughly into applications of the 
existing consistency test to more general networks and 
signals, and extensions and improvements to the algo- 
rithm. 

We will systematically study a larger variety of wave- 
forms than the small set considered in the paper. These 
should include two-polarization GWBs. (The supernova 
waveforms used here are linearly polarized). This wider 
set may include various supernova catalogs [8, 9, 10], and 
approximate waveforms for black-hole binary coalescence 
(see for example [7]). The latter are particularly impor- 
tant, since it is quite plausible that black- hole binaries 
will be the first transient signals to be detected. 

Another near-term goal is to apply our modified null- 



FIG. 9: Receiver- Operator Characteristic (ROC) plot for 
min(^nuii/^inc) as a statistic for distinguishing GWBs from 
noise glitches. The ROC curve value is given by the fraction 
of GWBs (true acceptance) and glitches (false acceptance) 
of given rms SNR falling to the left of a line of constant 
^nuii/^inc (a diagonal line) in Fig. 6. The rapid rise of the 
curves at low false acceptance is indicative of the ability of 
the method to confidently distinguish a significant portion of 
GWB signals from the glitch population. 



stream test to other networks, such as those with four 
and five detectors. A fourth non-aligned detector will 
reduce the fraction of sky over which only two detectors 
have significant sensitivity [21], increasing the strength of 
the test. An additional aligned detector, such as the two- 
kilometer detector at LIGO-Hanford, would also provide 
a second null stream without extra sky coverage. Laz- 
zarini et al. [ ] have demonstrated how the output of the 
two LIGO-Hanford detectors can be combined to form a 
single pseudo-detector with greater sensitivity than ei- 
ther. The difference between the detector outputs is also 
a null stream, effectively a detector with zero antenna 
response, which can be employed in our consistency test. 
This is a computationally cheap test since the H1-H2 null 
stream is independent of the sky position, providing the 
basis of a simple hierarchical analysis scheme. 

We also plan to test the power of our consistency test 
on real data, with real noise transients. We should note 
that the artificial noise coincidences studied in this work 
are pathological in the sense that they are injected with 
time delays and amplitude responses consistent with ac- 
tual sky positions, and waveforms that are identical (in 
individual detectors) to GWBs. Although we would not 
expect noise in actual detectors to be so pathological, it 
is important to characterize the robustness of our tech- 
nique for real data, and this was the motivation for our 
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FIG. 10: Receiver- Operator Characteristic plot (ROC) for 
min(£^nuii — ^inc) as a statistic for distinguishing GWBs from 
noise ghtches. The ROC curve value is given by the fraction 
of GWBs (true acceptance) and glitches (false acceptance) 
of given rms SNR falling to the left of a line of constant 
^nuii — ^inc in Fig. 7. The rapid rise of the curves at low 
false acceptance is indicative of the ability of the method to 
confidently distinguish a significant portion of GWB signals 
from the glitch population. This is not as powerful a statistic 
as min(^nuii/^inc) (see Figure 9). 



implementation. 

Another important aspect of real data is that their 
calibration could be inaccurate, implying an imperfect 
cancellation of the GWB by the null stream. We will 
study quantitatively the effects of realistic calibration er- 
rors on the effectiveness of the null-stream technique (see 
also the study by Ajith et al. [ ; ]). 

A first improvement is to optimize the algorithm for 
implementing our test. An optimal identification of the 
integration time and frequency band over which the null 
and incoherent energies are calculated will improve the 
effectiveness of the test by minimizing the amount of 
noise included. 

An important augmentation of the coherent analysis 
is improved techniques for determining the sky location 
of the source. Current efforts estimate the sky position 
as the extremum in a sky map of the null stream or 
some likelihood statistic [13, 20]. For example, Giisel 
and Tinto [ ] demonstrated that locating the minimum 
of the null energy allows one to determine the direction 
to the source of a kiloHertz GWB with high accuracy for 
typical SNRs of ~ 40 — 60 (converted to our Prms)- How- 
ever, as we have seen, sky maps exhibit structure which 
is a combination of the network geometry and the signal 
waveform. For example, linearly polarized signals pro- 



FIG. 11: Receiver- Operator Characteristic (ROC) plot for 
min(£^nuii) as a statistic for distinguishing GWBs from noise 
glitches. The ROC curve value is given by the fraction of 
GWBs (true acceptance) and glitches (false acceptance) of 
given rms SNR falling to the left of a line of constant ^nuii 
(a vertical line) in Fig. 8. The lower slope of the curves at 
low acceptance is due to the fact that many glitches have sky 
positions for which the null energy is small, even when the 
SNR of the glitch is high. The null energy alone is therefore 
not very effective for confidently distinguishing GWBs from 
noise glitches. While one could chose a threshold that varies 
with Einc to get better performance, min(£^nuii) is still not as 
powerful a statistic as min(^nuii/^inc)- 

duce interference fringes in the sky map. These fringes 
appear as rings of fixed time delay with respect to the 
various detector baselines. A global analysis which takes 
account of this structure could in principle average over 
local noise fluctuations in the sky map to achieve an im- 
proved pointing accuracy for weaker signals. 

We will further investigate Bayesian interpretations 
and formulations of the null-stream technique, and com- 
pare the effectiveness of this approach against the pro- 
cedure presented in this work. We will also investi- 
gate the possibility of incorporating additional tests such 
as distribution-free (non-parametric) correlation tests. 
These could prove valuable when analyzing real data 
since they enforce known statistics even when the noise in 
the data follows an arbitrary noise distribution function. 



APPENDIX A: NULL STREAM PROJECTION 
OPERATOR 

In this section we derive explicit expressions for the 
projection operator which acts on the network data vec- 
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tor to produce null streams. This projection operator 
projects the data orthogonally to and . 

Let m and n be any orthonormal pair of vectors that 
span the space spanned by and F^ (the vectors 
and F^ are not necessarily orthogonal). Then the pro- 
jection operator that removes the gravitational- wave con- 
tribution is 



D 



/3=1 



(Al) 



For example, choosing 



m 



F+ 

w 
I w I 

F^-{m-F^)m 
\F^-{m-F^)m\ 



one finds 



:>NS 



1 



WQ! WQ! 



-M. 



+ x 



det(M) 

where the 2x2 matrix M is defined by 

M = Fw^Fw . 
In matrix notation (A4) becomes 

F F ^ 



Mx+ M++ 



pNS ^ J 



(A2) 
(A3) 



F+ 
(A4) 



(A5) 



(A6) 



If F^ oc F^ (e.g., for co-aligned detectors) then the 
projection operator simplifies to 



dNS 

a/3 



{Sap - rriamp) 



(A7) 



with rria given by (A2). 



APPENDIX B: MINIMUM- VARIANCE 
WAVEFORM RECONSTRUCTION 

In this section we derive explicit expressions for the 
two amplitude components of the wave, /i+ and h 
optimally reconstructed as linear combinations of the de- 
tector data. This is a generalization of the technique used 
by Giirsel and Tinto [13], and was first derived by Flana- 
gan and Hughes [6]. 

We will assume that the whitened data streams 
have been time-shifted before these combinations 
are constructed, and will simply write d^aifi] 

Our goal can be reduced to the problem of identifying 
two vectors V+, V x such that 

hf = V+.dw = /i+ + V+.nw (Bl) 
/if = Vx • c^w = /^x + Vx • nw . (B2) 



From the above equations we conclude that these two 
vectors must satisfy the constraints 

V+-F^ = 1 

V+ ■ =0 

Vx •F+ =0 

V^-F^ = 1. (B3) 

This implies that the errors in estimating the two wave- 
forms are 

5h^f = h^^ -h^ = V^'h^ (B4) 
J/if = /if - /ix = Vx • nw . (B5) 

If the noises affecting the detectors are Gaussian dis- 
tributed, the "optimal" choice for V+,Vx is that im- 
plied by minimizing the mean-square errors of the noises 
affecting the two reconstructed waveforms. 



N-l 



k=0 
N-l 



D 



(B6) 



(B7) 



fe=0 



That is, we want the shortest (in the Cartesian sense) 
vectors V+, V x which satisfy the constraints (B3). This 
minimization can be performed by using the method of 
the Lagrange multipliers. For instance, for V+ we have 



= 



w/3 



(B8) 



where (A+,Ax) are the two Lagrange multipliers. The 
above equation leads to 



1 



(B9) 



Imposing the constraints (B3) in the above equation re- 
sults in the linear system of equations 



(BIO) 



where the symmetric 2x2 matrix M is defined in (A5). 
Inverting it follows that 





-Mx 



det(M) I M+> 



which gives 
V+ = 



det(M) 



(Bll) 



(B12) 
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By performing an analogous calculation for we find 
1 



det(M) 



(-F+,M+x+f„"«M++) . (B13) 



Note that the weighting of a particular detector da van- 
ishes if = = Fwa, which occurs if detector da 
has no sensitivity to the sky location being considered 
(F+ = = F^)^ or if it is much noisier than the other 
detectors {Sa oo). In these cases, the expressions 
(B12,B13) at that frequency and sky position reduce to 
those for the network that does not include detector a. 

In matrix form, V^'^ are the rows of a 2 x matrix 
V where from (B12), (B13), and (A5) 



V = M-^F^ = {F^F)-^F^ 



(B14) 



That is, the waveform reconstruction matrix V is the 
pseudo-inverse of F [ ]. Indeed, it is easily verified that 



VF = 1. 



(B15) 



This is equivalent to the conditions (B3). 

If F^ (X F^ (e.g., for co-aligned detectors) then only 
a single linear combination of the two polarizations h-^^x 
can be extracted. Choosing the polarization gauge t/j' 
such that K+l = + |F^^|2)i/2, ir;x = q, the 

reconstruction matrix reduces to simply 



V = 



W 

F'+P 



(B16) 



See Rakhmanov [31] for a discussion of singularities in 
waveform reconstruction due to rank-deficiency of F^, 
and the use of Tikhonov regularization to avoid such 
problems. Rakhmanov's expressions correspond to ours 
with the replacement F^'^ F^^^; i.e., replacing 
the antenna responses by the noise- weighted antenna re- 
sponses. 
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